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In this letter, we consider the constraints imposed by the recent measurement of 
B s — B s mixing on the new physics contribution to the rare decay B s — > fi + fj,~ . New 
physics in the form vector and axial-vector couplings is already severely constrained 
by the data on B — > (K, K*)[i + fi~ . Here, we show that B s — B s mixing data, together 
with the data on K° —K° mixing and Kl — > fJ- + ^~ decay rate, strongly constrain the 
scalar-pseudoscalar contribution to B s — > /i + \i~ . We conclude that new physics can 
at best lead to a factor of 2 increase in the branching ratio of B s — > fj, + [i~ compared 
to its Standard Model expectation. 

The flavour changing neutral interaction (FCNI) b — > sfi + fi~ serves as an important probe 
to test the Standard Model (SM) and its possible extensions. This four fermion interaction 
gives rise to semideptonic decays B — > (K, K*)n + n~ and also the purely leptonic decay B s — > 
u + u~. The semideptonic decays B — > (K,K*)u + u~ have been observed experimentally 

fifin, nnn, 

111, 12, 13| with branching ratios close to their SM predictions [4|, |5|, |6fl. At present there is only 



an upper limit, 1.0 x 10~ 7 at 95% C.L., on the branching ratio of the decay B s — > fx + fi~ [3, [8|. 

n 

The SM prediction for this branching ratio is (3.2 ±1.5) x 1CT 9 [9] or < 7.7 x 10~ 9 at 3a level. 
B s — > fi + \T will be one of the important rare B decays to be studied by the experiments 
at the upcoming Large Hadron Collider (LHC). We expect that the present upper limit will 
be reduced significantly in these experiments. A non-zero value of this branching ratio is 



lflL 



measurable, if it is > 10 

In a previous publication [ill], we studied the constraints on new physics contribution to 
the branching ratio of B s — > coming from the experimentally measured values of the 

branching ratios of B — > (K, K*)fi + fi~ . We found that if the new physics interactions are in 
the form of vector/axial- vector operators, then the present data on B(B — > (K, K*)fi + fi~) 
does not allow a large boost in B(B S — > n + n~). By large boost we mean an enhancement 
of at least an order of magnitude in comparison to the SM prediction. However, if the new 
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physics interactions are in the form of the scalar/pseudoscalar operators, then the presently 
measured rates of B — > (K, K*)fi + fi~ do not put any useful constraints on B s — > /x + /x~ and 
Bnp(B s —> /x + /x~) can be as high as the present experimental upper limit. Therefore we are 
led to the conclusion that if future experiments measure B s — > /x + /x~ with a branching ratio 
greater than 10~ 8 , then the new physics giving rise to this decay has to be in the form of 
scalar/pseudoscalar interaction. 



Recently B s — B s mixing has been observed experimentally [13| . with a very small exper- 
imental error. In this paper, we want to see what constraint this measurement imposes on 
the new physics contribution to the branching ratio of B s — > ji + ji^ . In particular, we con- 
sider the question: Does it allow new physics in the form of scalar/pseudoscalar interaction 
to give a large boost in Bnp(B s — > /i + /x~) ? 

We start by considering the B s — > /i + /x~ decay. The effective new physics lagrangian for 
the quark level transition b — > due to scalar/pseudoscalar interactions can arise from 

tree and/or electroweak penguin and/or box diagrams. We parametrize it as 

= Gi Kaf + 9pi 5 )s fibfir + g^h, (1) 

where G% is a dimensional factor characterizing the overall scale of new physics, with di- 
mension {mass)~ 2 . This factor essentially arises due to the scalar propagator in tree or 
electroweak penguin diagrams (or scalar propagators in box diagrams) which couples the 
quark bilinear to the lepton bilinear. gf P and g^p are dimensionless numbers, charac- 
terizing, respectively, b — s and /x — /x couplings due to new physics scalar/pseudoscalar 
interactions. Electromagnetic penguins necessarily have vector couplings in the lepton bi- 
linear so they do not contribute to the effective lagrangian in eq. jl]). The amplitude for the 
decay B s — > l + l~ is given by 

M(B S - /x+/x-) = G, g s P b (0 \b l5 s\ B s ) [g^u(p^)v( Pfi ) + g^uip^vipp)} . (2) 

The pseudoscalar matrix element is, 

(0 b l5 s B.) = -i^^, (3) 

where m b and m s are the masses of bottom and strange quark respectively. 
The calculation of the decay rate gives 
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We see that the decay rate depends upon the new physics couplings (g s p) 2 and GfKg^) 2 + 
(gp 1 ) 2 ]- To obtain information on these parameters, we look at B s — B s mixing together 
with Kl — > decay and K° — K° mixing. 

Let us consider B s — B s mixing to obtain a constraint on {gp b ) 2 . Replacing leptonic 
bilinear by quark bilinear in eq. [TJ we get AB = 2 Lagrangian, 

l S bIb s = G 2 b(gf + g s P b ^)s b(9s + 9 s h 5 )s, (5) 

where G 2 is another dimensional factor. As in the case of G\, introduced in eq. (DQ), G 2 also 
arises due to the scalar propagator (or progators in the case of box diagrams) . Therefore it 
also has dimension {mass)~ 2 and is of the same order of magnitude as G\. From eq. ©, 
we calculate the mass difference of the B s mesons to be 

A% s = J^G 2 (g$) 2 B B . (6) 



2M Bs "'(m 6 +m a ) 2 
Thus the effective b — s pseudoscalar coupling is obtained to be 

b a = Am B3 (m b + m s ) 2 
[gp> 2B Bs f B M B G 2 ■ 1 ' 

We now consider the decay K L — > /i + \i~ . The same new physics leading to the effective 
b — > s[i + [i~ lagrangian in eq. (Dp), also leads a similar effective lagrangian for s — > dfi + fi~ 
transition. The only difference will be the effective scalar/pseudoscalar couplings in the 
quark bilinear. Thus we have, 

LtZ^r- = Gi s(9s d + gfi^KgT + nF-tii*. (8) 

The calculation of decay rate gives 

T NP (K L - mV) = 2(^) 2 [(^) 2 + ^g ( J + ; )2 - (9) 

Here extra factor of 2 occurs because the amplitudes A(K° — > /i + /i~) = A(K° — > fi + fi~) 
and K L = We see that G?[(^) 2 + (g 1 ^) 2 } can be calculated from T(K L -> 

once we know the value of (g s p) 2 - In order to determine the value of (g s p) 2 , we consider 
K° — K° mixing. The effective scalar/pseudoscalar new physics lagrangian for this process 
can be obtained from that of s — > dfi + /j,~ by replacing lepton current by corresponding quark 
current or equaivalently from effective lagrangian of eq. ((5j) where b — s quark bilinear is 
replaced by s — d quark bilinear, 

= G 2 s(g s s d + g s P d l5 )d s(gf + gf l5 )d. (10) 
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From this lagrangian, we obtain the Kl — Ks mass difference to be 



-G 2 (g P s ) B K 



PkM* 



2M K "(m s + m d y 

Thus the effective s — d pseudoscalar coupling is 

sd 2 2Am K (m d + m s f 



(11) 



Substituting the above value of (g s P d ) 2 in eq 
G\[{g s 



B K P K M\G 2 
I, we get 



(12) 



2r, ../'/' ^ + Ui'l',' f — 2 7r G 2 B K 



M 2 K Am K 



T NP (K L -> fx + /ji~ 



(13) 



Substituting the value of Gj[(g^) 2 + (g^) 2 } from eq. ([T3j) and (gp b ) 2 from eq. ((3) in eq. 

we get 



NP 



{B s 



fJL fJL 



M, 



K 



( Ami 



Am 



K 



B 



K 



B f 



NP 



fX fJL 



The branching ratio is given by, 



Bnp(B s 



1 (M B A 2 /Am B , 



B 



K 



r(B s ) 
t{K l ) 



B NP (K L -> 



(14) 



(15) 



2 \M K J \Am K J \B Bs 

We wish to obtain the largest possible value for B(B S — > ji + //~). To this end, we make the 
liberal assumption that the experimental values of Auib s , Amx and Bnp(Kl — > /i~V~) are 
saturated by new physics couplings. The decay rate for Kl —> n + jjT consists of both long 
distance and short distance contributions. The new physics we consider here, contributes 
only to the short distance part of the decay rate. In ref [l^j], an upper limit on the short 
distance contribution to B{K L — > is calculated to be 2.5 x 10~ 9 . The mass difference 

of the B s mesons is recenly measured by the CDF collaboration to be Am Bs = (1.17±0.01) x 
10~ n GeV [l^|. The bag parameters for the K and the B s mesons are Bk = (0.58±0.04) and 
Bb s = (1.30±0.10) [15l|. The values of the other parameters of eq. ( TT5i ) are taken from Review 



16|]: Am K = (3.48±0.01)xKT 15 GeVT(£ s ) = (1.47±0.06)xlO- 12 Sec 



of Particle Properties 

and t(K l ) = (5.11 ±0.02) x 10~ 8 Sec. Substituting these values in eq. ([TBI) , we get 



B NP {B S -> = (6.30 ± 0.75) x 10" 



where all the errors are added in quadrature. At 3cr, B S m(B s 



(16) 



< 7.7 x 10~ 9 where 



as Bnp(B s — > ) < 8.55 x 10 9 . Thus we see that this upper bound is almost the same 
as the SM prediction even if we maximize the new physics couplings by assuming that they 
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saturate the experimental values. Therefore the present data on B s — B s mixing together 
with data on K° — K° mixing and Ki — > decay puts a strong constraint on new 

physics scalar/pseudoscalar couplings and doesn't allow a large boost in the branching ratio 
of B s — > fi + fi~. 

We now assume that the new physics involving scalar/pseudoscalar couplings accounts 
for the difference between the experimental values and the SM predictions of Am^, Am Ba 



and 



the s 



by HQ 



(Am Ba 



lort distance contribution to T(Kl — > yU + /i ). The SM value for B s — B s is given 



>SM 



6tc 2 



Vb M Bs (B B Jl)M^S{x t )\V ts \ 



[1.16 ± 0.32) x 1(T U GeV, (17) 



with f Bs y B Bs = (262 ±35) Me V [15fl , rj B = 0.55 ±0 



).oi0 



and \ Vt 

n 



ts 



0.0409 ±0.0009 



S(x t ) with x t = rn\hmjyy is one of the Inami-Lim functions |1S|. The SM value for K° - K° 



mixing is given by 
(Am K ) SM 



20fl 



Or F 

6^2 



B K f 2 K ) M K M^ \* c %S(x c ) + X* t %S(x t ) + 2\* c \* tV3 S(x c , x t 



(18) 
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where Xj = V* s Vjd, Xj = m 2 Jm'^ ¥ . The functions S are given by 

c s r, x t 3x t 3xtlnx t 
S{x c ,x t ) = x r In — 



(19) 



x c 4(1 -XA 4(1 -x t ) 



Using = (L32 ± 0.32) [23j, r/ 2 = (0.57 ± 0.01) [18J, m = (0-47 ± 0.05) [2J, |25fl, B K 



(20) 



(0.58±0.04) tl5fl ; f K = (159.8±1.5) MeV, \V CS \ = 0.957±0.017±0.093, \V cd \ = 0.230±0.011, 
\V U \ = 0.0409 ± 0.0009 and \V td \ = 0.0074 ± 0.0008 H, we get 

{Am K ) SM = (1-87 ± 0.49) x 10~ 15 GeV. (21) 

All the masses were taken from Considering only the short-distance effects, the SM 

branching ratio for K L — > fi + fi~ in next-to-next-to-leading order of QCD is (0.79 ± 0.12) x 
10- 9 H. Substracting out the SM contribution from the experimental values of Am Bs , 
ArriK and Bnp(Kl — > , we get 



Bnp{B s 



fJL fJL 



B K 



{Am Bs ) exp - {Am Ba )sM 1 ( M Bg 
_ (Am K )ex P - (Am K ) S M \ 2 V M K 

[B exp (K L -> shor t - B SM (K L -> fi + fi~) 



B 



Us 



r{B s ) 
r{K L ) 



(22) 
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Substituting the experimental values and the SM predictions in the above equation, and 
adding all the errors in quadrature, we get 

B NP (B S -> n+fT) = (0.08 ± 2.54) x 10~ 9 . (23) 

which is consistent with zero. At 3a, the upper limit on the new physics contribution is close 
to SM prediction. Thus the present data on ArriB s along with ArriK and B^p(Kl — > /i + yT) 
puts strong constraints on new physics scalar/pseudoscalar couplings and doesn't allow 
a large enhancement in the branching ratio of Bnp(B s — > fi + y~) much beyond the SM 
predictions. New physics at most can cause a factor of two enhancement but 
not an order of magnitude. Hence the total branching ratio which is the sum of SM 
contribution and new physics contribution will be of the order of 10~ 8 and hence reachable 
at LHC. 

Conclusions: 

In this letter, we considered the constraints on the New Physics couplings of 
scalar/pseudoscalar type in the b — > s transition. It was shown previously that only such 
New Physics can give rise to an order of magnitude enhancement of the decay rate for 
B s — > Using the recent data on B s — B s mixing, together with the data on K° — K° 

mixing and the short distance contribution to Kl — > we obtained very strong bounds 

on B(B S — ► yU + /i _ ). New Physics in the form of scalar/pseudoscalar couplings can at most 
increase the B(B S — > by a factor of 2 compared to its Standard Model prediction. 

An order magnitude enhancement, previously allowed, is ruled out. 
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